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Abstract 



finite-difference equation which is solved in WKB approximation for large conformal spins of the 



We identify an integrable one-dimensional inhomogeneous three-site open spin chain which 
^ ; arises in the problem of diagonalization of twist-three quark-gluon evolution equations in QCD 
in the chiral-odd sector. Making use of the existence of a non-trivial 'hidden' integral of motion 
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ON • the problem of diagonalization of the evolution kernels is reduced to the study of a second order 
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Q-h! three-particle system. The energies (alias anomalous dimensions) of eigenstates with different 
scale dependence are found in limiting cases and compared with numerical calculations. 
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1 Brief phenomenological profile and outline. 



Deep inelastic scattering of leptons off nucleons, £N — > £'X, served for a long time as the only 
and the most reliable piece of knowledge about the partonic structure of hadrons in high energy 
collisions via measurements of non- and polarized chiral-even structure functions [![]. However 
the feasibility of measurements of other inclusive characteristics of hadrons was strongly limited. 
Nowadays in order to study in detail the high-energy dynamics of constituents, especially in 
polarized targets, one is forced to use complementary sources of experimental information. In 
order to circumvent the above disadvantages of conventional deep inelastic scattering process 
one has to resort to nucleon-nucleon collisions, like the Drell-Yan lepton pair production hh' — » 
(7*, W, Z) — > £ + £~X, which are not affected by the above restraints and allow to explore a wealth 
of hadronic properties, e.g. to access a new type of parton densities — chiral odd distributions 
hi at twist-2, Hl and e at twist-3, and, related by QCD equation of motion to the latter two, 
generalized three-parton correlators 0. The twist-3 distributions being a direct manifestation of 
the quark-gluon interaction in non-perturbative domain could serve as a testing ground for the 
phenomenological models of confinement and, perhaps, future lattice simulations. 

Similar to the more familiar chiral-even transverse spin structure function g 2 , although being 
of higher twist turns out to be unique by its feature to appear as a leading unsuppressed effect 
in longitudinal-transverse asymmetry in p\\p± — > £ + £~ X reaction 0, [|. However, an inevitable 
complication of all experiments is that the data is taken at different values of a hard momentum 
transfer Q and, therefore, to deduce the shape of distributions or to check sum rules one needs 
to know the way to describe the variation of quantities in question when going up(down)wards 
with Q from the point of view of the underlying microscopic theory of strong interaction - 
Quantum Chromodynamics. Within the latter framework this change is described by evolution 
equations which in turn ctnsG £is du result of intrinsic ultraviolet divergencies of the product of local 
elementary fields, — which define an operator content of distributions, — separated by light cone 
distances, and thus require renormalization which introduces a scale dependence. The specifics of 
twist-3 sector as compared to the well-studied twist-2 case is that the above mentioned structure 
functions being defined as two-particle hadron-hadron matrix elements receive contributions from 
the three-parton quark-gluon correlators 0, [|]. Resolving this complication in favour of a study 
of the renormalization group evolution of the latter one ends up with a very complicated problem 
of working out the mixing of three-particle local operators which in general cannot be handled 
analytically. 

One can ague that it is not really possible to see this logarithmic violation of scaling exper- 
imentally due to the contamination of a more strong 1/Q background behaviour of prefactors 
which accompany the asymmetries. Nevertheless this knowledge is important for contrasting the 
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theoretical models with the measurable quantities. Since the model estimations are usually done 
at a rather low momentum scale the evolution effects are very prominent and change in a siz- 
able way the shape of an input distribution during short evolution times. So that perturbative 
evolution can be thought of as a prime mover of gross features of non-perturbative functions at 
hadronic scales, however much larger then Aqcd- Recall in this respect the radiatively generated 
parton densities 0. 

The distribution functions Hl and e as measured in most of processes are given, up to 'kine- 
matical' Wandzura-Wilczek-type contribution in Hl which obeys the ordinary twist- 2 evolution, 
by particular integrals (see Eqs. (|2],|3|) below) of three-particle quark-gluon correlation functions. 
The first drastic simplification of their evolution occurs for large number of colours. In multi- 
colour limit the weight in the integrand of convolution integrals coincides with the eigenfunction 
of the leading order evolution kernel and, therefore, the above mentioned chiral-odd structure 
functions evolve multiplicatively — the so-called 'decoupling' phenomenon [7|, |[ [|. However, in 
the spin-independent Drell-Yan process the three-particle correlators associated with the unpolar- 
ized distribution e enter folded with different coefficient functions || and, therefore, will not evolve 
autonomously even in multicolour limit. In the realistic QCD case, when the correction are 
accounted for, the structure of evolution becomes more complicated: all three-particle operators 
are of relevance. But the violation of the simple pattern of multicolour evolution can be well 
understood once it is the case for large- N c limit. Since these corrections produce oc (10 — 20)% 
effects in the evolution of the model distributions^ which although being small and apparently 
might be sufficient for practical purposes at present should still be under control. This paper is 
devoted to an analysis of the large- iV c spectrum of the three-particle evolution equations. 

Our consequent presentation is organized as follows: In the next section we describe an ap- 
proach to study the scale violation of the so-called quasi-partonic operators which are related to 
the quantities in question by QCD equations of motion. We show that the tree-level conformal 
symmetry of classical QCD Lagrangian allows for a partial diagonalization of the evolution ker- 
nels. In section 3 we elaborate on this point further and construct an appropriate three-particle 
conformal basis. A distinguishing property of the problem at hand is its equivalence to an exactly 
solvable one- dimensional open spin chain which is established in section 4. Making use of the 
existence of a 'hidden' integral of motion allows to reduce the eigenvalue problem for the energy of 
the system (alias anomalous dimensions) to the solution of a second order finite-difference equa- 
tion — referred to as the master equation — which is derived in section 5. Consequent discussion 
concerns the study of the analytical structure of the spectrum of the conserved charge for large 

2 We discard from consideration some pathological shapes of gluon distribution in three-parton correlators which 
induce non-smooth evolution ItLll. 
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values of the total conformal spin. Its asymptotic solutions are given in section 6. These results 
allow to describe fairly well the spectrum of energy eigenvalues found by numerical diagonalization 
of the mixing matrix of the anomalous dimensions. Finally, we go to conclusions. 



2 Approach to scaling violation phenomena. 

In QCD the twist-3 distributions, Hl and e, are expressed via the light-cone Fourier transformation 
of bilocal quark string operators || 

{£!)) = f /^W(°)*H (JJ *MI*>. « 

with a phase factor $[0, nn] which ensures the non-abelian gauge invariance. Although these defi- 
nitions resemble the ones of ordinary twist- 2 parton densities a more closer inspection immediately 
reveals an essential difference from the latter: twist-3 functions involve the "bad" light-cone com- 



ponents of the fields in the language of the Kogut-Soper infinite momentum frame formalism [[TO 
By virtue of the QCD Heisenberg equation of motion for the quark field and constraints coming 
from Lorentz invariance, the functions (|T]) can be rewritten in an explicitly interaction dependent 
way asQ 

e(x) = J^Z(x,P), (2) 



hi{x) = 2x 2 ' 



r 1 dp 


f dp' 


' d 


d ' 


L ~w j 




dp 





Z(P,{3'). (3) 



!x p 2 

in terms of correlation functions 

Z( Xl , x 2 = - Xl - x 3 , x 3 ) = I ^L^*e***i*«**(h£ 0, ks) ± Z(-«3, 0, -kj)} \h) (4) 

of the non-local quark-gluon operators 

«i,«2,«s) = ^(«3n)oi. [ 1 J t a gG"^(K 2 n)^(K 1 n), (5) 
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with gauge link factors omitted for brevity. However the operators in Eq. ([5]), while contain one 
parton more, share a bulk of properties with the usual twist-2 operators: only "good" light-cone 
components and collinear momenta appear. "Bad" fields and transverse momentum dependence 
have been eliminated in favour of interaction which originates them. Thus in this description 
the theory is endowed again with a parton-like interpretation albeit of the QCD improved quark 



3 The first contribution on the r.h.s. of Eq. (H) is afore mentioned 'kinematical' twist-2 WW-term 
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Figure 1: The diagrammatical representation of the generic form of the leading order evolution 
equation for the three-parton correlator Z(xi,X2,x 3 ) in the light-cone gauge. The blob K stands 
for the pair-wise kernel of interaction of t-channel particles. 



model. And these are the properties which distinguish a special class of the so-called quasi-partonic 
operators (T^j which are the most appropriate for the study of renormalization group dependence 



as the corresponding evolution equation turns out to be homogeneous. This latter property has 
arisen due to the fact that the pairs of parton fields in Eq. ([5]) correspond to the familiar twist- 2 
operators although they can have fermionic quantum numbers besides bosonic. 

2.1 Evolution equation. 

The generic form of the evolution equation for twist-3 correlation functions in the momentum 
fraction formalism looks like (see Fig. [I]) 

(J> 2 -^Z(X 1 ,X 2 ,X 3 )= J Il dx i S i^ X '^j K ({ Xi }\i X 'i}) Z ( X 'l> X 2> X '3)> ( 6 ) 

with the constraint x\ + x 2 + x 3 = of the forward inclusive kinematics for the momentum 
fractions of partons imposed on the both sides of this equation. Due to topology of one-loop 
Feynman diagrams the evolution kernel K has a simple pair-wise structure at leading order: 

K = J2 K ( x i, x j\ x 'n x 'j) $ (xi + xj -x[ -x^) , (7) 

i<j 

with K '■ ) being a twist-2 interaction kernel of two nearby particles. 

The solution to the leading order evolution equation (|Bp can be written in the form 

Z(x u x 2 ,x 3 \Q 2 ) = ^{a } (xuX2,x 3 )(^^) ((Z {a M))) (8) 

{a} \°isW )/ 

of expansion w.r.t. polynomials \Ps a \{xi, x%, x 3 ) which are the eigenfunctions of the leading order 
kernel (|7|) with eigenvalues £{ a y and where ((Zi a y)) denote reduced hadronic matrix elements 
of local operators defined as (Mellin-like) moments of the three-parton correlation function. The 
summation goes over a set of quantum numbers {a} which enumerate the solutions, e.g. spin of the 
system etc. And /3 = ^T F Nj — y Ca is the first coefficients of the QCD Gell-Mann-Low function. 
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Since the polynomials depending on the momentum fractions can be immediately translated by 
Fourier transform to the ones depending on derivatives acting on the elementary field operators, 
they will correspond to the eigenvectors of the anomalous dimension matrix in a basis of local 
operators. Thus the problem of solving the evolution equation is equivalent to the diagonalization 
of the matrix of anomalous dimensions 

where we have explicitly extracted the dependence on the coupling constant a s . 

For the problem at hand the total evolution kernel takes the same form for the parity-even (e) 
and -odd (h L ) sectors 

K {{xiWix'j}) = ~^{ qqK (8) (x 1 ,x 3 \x' 1 ,x' 3 ) 5(x! +x 3 -x[ 

+ q9 K(3) (xi, x 2 \x' 1 , x' 2 ) 8{xx + x 2 - x[ 
+"K^ (x 3 , x 2 \x' 3 , x' 2 ) 5(x 3 + x 2 - x' 3 

with known non- forward evolution kernels, ^ 1( ^ 2 K^, for colour triplet (c = 3) vector (D = V) 
and colour octet (c = 8) transversity (D = T) sectors. The selection of particular solutions 
comes from different parity properties of three-variable functions Z and Z w.r.t. permutation of 
quark fields: the former is symmetrical Z{x\,x 2 ) = Z(x 2 ,Xi) while the latter is antisymmetrical 
Z(xi,x 2 ) = —Z(x 2 ,xi) under this interchange. We have added in Eq. ( |T0"1) the trivial charge 
renormalization counterterm which is traced back to the presence of the coupling constant in the 
definition of the composite operator (|5|). 



x. 



(10) 



4) 



Ob Q 



^<5(xi - x\)5(x 2 



2.2 Eigenvalue problem. 

Due to absence of the QCD /5-function contribution in leading order pair-wise evolution kernels 
we can use the conformal invariance property of QCD at tree level since the renormalization 
counterterms of composite operators at leading order are just the ones in classical theory while 
from the side of QCD Lagrangian the only counterterms which are required in the physical sector 
to cure the divergencies are the field renormalization constants. But the anomalous dimensions 
associated to them can be embedded into redefined representations of the conformal group - 
replacing the canonical by the scale field dimensions — and thus preserve conformal covariance. 

Making use of this property of classical QCD one can diagonalize the pair-wise non-forward 
kernels |13|-||18| and present them in a factorized form |19| 

2 J^0 J \Xi+X 2 J J \X 1 + X 2 J 
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where w{x x ,Vx\x 2 ,v 2 ) = x^x? and uj j (p 1 ,p 2 ) = ^j+^+^+mru-^+li+i) ^ with v t = <b + st- 1, 
and dn and standing for the canonical dimension and spin of the constituent (f>t = {q,g}- The 



p{v*,vi) are usua j j aco bi polynomials p0| which form an infinite dimensional representation of 



the collinear conformal group pi| , — whose su(l, 1) algebra 

[.7 3 , J*]. = ± J*, JT]_ = -2J 3 , (12) 

is defined by generators of the momentum J + = iP+, special conformal transformation J~ = |/C_ 

and particular combination of dilatation and angular momentum J 3 = |(P + A4 (-), — in the 

space spanned by bilinear (in elementary field operators) local operators with total derivatives 

H-flU 



O jk = <p 2 {idi + iarfFf*** [^fj 4>i- (13) 

These operators possess conformal spin J\ 2 = j+^(vi+v 2 +2) which is the eigenvalue [J7 -2 , Ojk]- = 
J 12 (J 12 - l)O jk of the Casimir operator J 2 = J 3 (J 3 - 1) - J + J~. 

In the Eq. ( [TT| ) , all dynamics of the system is concentrated in the anomalous dimensions of the 



8,-d- 



conformal operators, ^^j. Thus the former depend only on the eigenvalues of the two-particle 
Casimir operators in subchannels. For the eigenvalue problem (§) it can be translated to the 
fact that the pair-wise kernels commute with generators of conformal transformation, thus, again 
being a function of Casimir operator only. 

This dependence can be easily traced from the corresponding eigenvalues. For our consequent 
considerations we need the bosonic colour-octet chiral-odd, q ^ljs) ■ anc ^ fermionic colour-triplet 
chiral-evenF] 95 7(3) . anomalous dimensions which read 

q %), = (C F - + 2) - #(1) - 3) , (14) 



2 

7(3) j = ^(2^(j + 1) + 2^(j+4)-4^(1)) 



+ ( c '-t) (, + A + ,) + ^-^-M 

where a(j) = (— 1) J . 

Therefore, the evolution equation can be reformulated into the eigenvalue problem for the 
three-particle system in a basis of local operators 



H qcd V = £ Q cd&, (16) 



where K = —%^7i with Hamiltonian 



n QCD = ^f 8) (Ji 3 ) + ^) + q ^U-^) ~ 



We do not bother about parity since it is irrelevant for fermionic operators. 
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and we have defined qq jj^2 = 2 qq H(J) and qa jj-5/2 = 2 qg H(J). The operators J's are formally 
determined as solutions of the operator equationF] for Casimir operator in some representation^] 
f = J(J-l). 

One of the complications to find the solution of the eigenvalue problem (|16j) is a rather non- 
trivial dependence on the colour structure of the quark-gluon system. In order to get rid of it we 
take large N c limit in Eq. (|17"D which results into omission of the l/iV c 2 terms in Eq. flTJ). This 
leads to the Hamiltonian 

H= Y { kl2 + ~ 5} ' ^ h ab = ^(j<*-^+$(jab+^-2ll)(l), (18) 

whose diagonalization will be our cherished goal. 



3 Conformal basis: General formalism. 

As we have seen above the conformal covariance provides a partial diagonalization of the three- 
parton evolution kernels. Therefore, an appropriate basis constructed with regards to the tree-level 
conformal symmetry can simplify the problem further. This question is considered below in the 
present section. 



3.1 Auxiliary #-space. 

It turns out to be convenient to introduce an auxiliary vector space where the conformal invariance 
properties are reformulated in a manifest and simple way. This description appears once we 
consider the elements 

e k = (19) 

as building blocks of the formalism. In this representation [J' ± ' 3 , x{&)]- — J ±,3 x{@) t ne generators 
take the following form 

J+ = ( v + iy> + fi±, J- = ±, j.= i(„ + i ) + s », (20) 

and the quadratic Casimir operator is given by 

f = J3(j3_!)_ j+j-. (21) 



°Note that the operators J introduced here differs from the ones defined in Ref. |1 by a shift: J horc = + 1. 
6 Once conformal invariance is established we can abstract from the particular representation used for the 
derivation. 
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For a multi- variable function x(9i, &2, ■ ■ ■ , 9 n ) the operators are defined by the sum of single particle 
ones as J ±,z = J2"=i Jp 3 an d they obey the usual commutation relations: [J 3 , J^- = iJ*, 
[J + , J~\- = — 2J 3 . Obviously, the single particle state 0% is an eigenstate of the Casimir operator 

f e 9t = f t 9 e , with = (22) 

Since the Eq. (|16|) defines the physical anomalous dimensions, the latter have to be real that 
requires the Hamiltonian (|18|) to be selfadjoint[] w.r.t. an appropriate scalar product. We define 
the latter as 

( x (9 1 ,9 2 ,...,9 n )\x(9 1 ,9 2 ,...,9 n ))= f dM X (9~i, ■ ■ ■ , 9 n )x(0i, fc, • • • , 9 n ), (23) 

J a 

where 

cM^n^a-W^ 1 and Q=\J{\9 e \<l} 

i=\ Zm 1=1 

and 9 = 9*. Then it can be seen that the Casimir operator and, as a consequence, the Hamiltonian 
are selfadjoint operators w.r.t. such defined inner product 

( j 2 )t = f, vO = n, 

and therefore Im£ = 0. 



3.2 Fundamental basis. 

Using the above definitions one can construct a second order differential operator (here and 
throughout 9 ab = 9 a — 9 b ) 

J\ = -t^Wj 1 ^ + \^ - (d a + d b ) , (24) 

which corresponds to the two-particle Casimir operator in (a, 6)-subchannel. Its eigenstates 
(jvz) 9i2 coincide with bilinear conformal operators (13) and possess the same eigenvalues. The 
three-particle Casimir is an obvious generalization 

3 

s2 A 2 A 2 a 2 ^ — ^ ^2 / \ 

J = ^ 12 + J 23 + *^13 — 3t ■ (25) 

i=i 

As it has been advocated in Refs. P^|-[pi]| it proves convenient to introduce the following 

basis of three-particle operators which simultaneously diagonalize the total three-particle Casimir 

7 To be precise Hermitean matrices do possess real eigenvalues, however, the inverse statement in not true. Recall 
in this respect triangular matrices with real diagonal elements. 
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-2 "2 

operator J and one in a subchannel J abJ say a = 1 and 6 = 2, 

S 2„ ,„ „ ,„ x / , 1 , A / , 1 



?VjyiPuh%) = (^+2^1 + ^ + ^3 + 3)) (j + -(^l + ^ + ^ + l))Pj;i(^1^2|^), (26) 

J^J^iAl^) = (j + 2>i + ^)) (j + ^ + ^+2))^(0 1 ,0 2 |0 3 ). (27) 



The states with unit norm w.r.t. the scalar product (p3|) are defined immediately from the above 
equations and read 



k 



Vj-AOi, W = M-\J,M, I*, E ^4 ^ fc | 9 , ^ fc 0%)X (28) 

^ (2j + v x + z/ 2 + 2) fe v ' 



with normalization factor 



JK-J - J')!(J + + ^2 + l)j {J + 3 + V1 + V2 + V3 + 



A/ ^2, ^3) = — 7 — ^ l—s -7—s 777— r—r, -, (29) 

{Vl)j+l{V2)j+l{V3)j-j+l (2j + V\ + v 2 + 2)j_ j 

where we have introduced the Pochhammer symbol (a)e = T(a + £)/T(a). By construction the 
above basis is orthonormal (Vj'-f(9i, 9 2 \9 3 )\Pj- 1 j{9\, 9 2 \9 3 )) = Sj'jSj'j- 

3.3 Matrix elements. 



The form of the basis given by Eq. (28) is especially suitable for evaluation of matrix elements 

of SU(1, 1) generators. The 'proper' generators of the basis V(9i, 9 2 \9 3 ) are diagonal (^,^) and 

-2 

the only operators which possess off-diagonal elements are J ^ 2 )3- The important property of the 
basis ( P8| ) which distinguishes it among other possibilities is that the matrix elements of 'improper' 
operators are three-diagonal only 

~2 / "2 -2 -2 \ 

( / P.j';j'(9 1 , 9 2 \9 3 )\J 23 \V J . j (9 1 , 9 2 \9 3 )) = 5j,j \ Sfj[J 23 ] ji: j + Sf ;j+1 [J 23 \ j+1 j + Sf,j-i[J 23 ]j-i,j ) (30) 
where the diagonal part is 

[Jlshi = \ {[f]jj ~ [fuh + E ft + til ~ fx) ([An - Si) [JI2UI) , (31) 

and the non-diagonal elements are 

rf 2 i Af(J,j>i,^ 2 ,^ 3 ) / ■ 1 1 w t .. v , QO x 

= MU 3 + l\v lM ) {3 + 1)( J - J + Us) - (32) 



From hermiticity it follows that [J 23 ]j+i,j = [J 



2 3 \j,j+i- 
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3.4 Other representations of the basis. 



Since the operator (p8|) is the highest weight of the conformal group in the three-particle basis, 
i.e. J'V = 0, and is the homogeneous polynomial (through Eq. (p6|)) of degree J by the Euler 
theorem, this suggests that the net dependence of Eq. ( |2"8| ) is on a single variable only (up to a 
prefactor). This is really the case and immediately one can perform the resummation of the series 
in Eq. fl28|) to obtain a compact expression for the #-space basis (cf. pc 



v J . j (9 1 ,e 2 \e 3 ) =Af-\j,j\v 1 ,v 2 ,v 3 )9( 2 9i- J 2 F 1 ( J 2 . + J ^^2 



9 ) , with 9 = (33) 

/ #32 



In general the solution to the condition J~V = is given by a translation invariant polynomial, 
e.g. V oc J2ji+j 2 =j c ji,j2^i2^32 with undetermined expansion coefficients. This is the standard choice 
which leads to Appel polynomials when transformed back to the space of local operators (TB| . If, 



as we have done above, we will impose an additional constraint (p7|) we get a recursion relation 
for c's which when solved gives (p3|). 

The relation to the space of local operators can be readily figure out by using the definition 
(IS) so that after some algebra we end up with (cf. 



r(j + ^ + i)r(j + ^ + i) 

r(j + i) 1 } 

v r( J - J + us + 1) r( J + j + y± + v 2 + 2) r ^ 

x r ( j- J + i) r(2 J + , 1 + , 2 + 2) ^Jh^^)^(gi.g»lg8)- 

3.5 Racah decomposition. 

Since the choice of the second condition ( p7|) on the basis is arbitrary, — one can analogously take 
a requirement for the element of the basis to be an eigenstate of operator J, x 2 \ 3 , — one poses the 
question as to how the transformation from one basis to another is accomplished. Obviously, we 
meet a standard quantum mechanical problem of addition of three 'angular momenta'. According 
to general theory they are related via the decomposition 

j 

Pv&Vj-MM0z)=Vj»{0zMh) = Y, w AJ)Vj; k (9u9 2 \9 z ) ) (35) 

fc=0 

where the coefficients Wjk(J) are viewed as Racah symbols for SU(1, 1). The operator P13 per- 
mutes the fields '1' and '3' in teh operators Vj-j. The explicit form of Wjk(J) can be easily 
deduced from the above formula sinceQ Wjk(J) = (Pj;k(9i, 9 2 \9 3 )\Vjj(9 3 , 9 2 \9i)) and making use 

8 The Wjk(J) are the matrix elements of the operator P13. 
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of the representation ( p4|) with the derivatives loosely interpreted as momentum fractions. Going 
to the special system with idi = z, id 2 = 1 — z and id 3 = 0, where the calculation methods of Ref. 



18] can be applied in a straightforward way, we end up with the result 

r(j 



Wjk(J) 



-i 



\ J — h 



T(J + k + u 1 + u 2 + u 3 + 2) 



x 



X 



X 



r( j - k + i)v(k - j + 1) r(j + j + v x + v 2 + v z + 2) 

r(fc + + l)T{k + v 2 + l)T(k + V1 + V2 + I) 
T(j + u 2 + l)T(2k + ui + u 2 + l)T(2k + v x + v 2 + 2) 

r( j + fc + z/i + v_2 + 2)r(2j + v 2 + 1/ 3 + 2) wy, fcj^, z/ 2 , ^ 

T(J -J+V1 + l)r(j + fc + l/ 2 + I/ 3 + 2) Af(J,j\v 3 , Z/ 2 , */l) 

„ / - J + k, J + k + v x + v 2 + v 3 + 2, k + v 2 + 1, k + 1 
j + k + v 2 + v 3 + 2, 2k + v x + v 2 + 2, - j + 1 



4-^3 



These coefficients are real = W* k and orthogonal 



(36) 



£=0 



They satisfy the three-term recursion relation 
^2 + ^3 



J + 



, 1/2 + 1/3 1 1 
J + — - — + 1 



W jk {J) 



23jfc,fc- 



(37) 



_ x W jik+1 {J) (38) 



with matrix elements of the Casimir operator «7 23 evaluated above. Using it one can easily derive 
useful approximate formula in the limits r = j = fixed, j = fixed, J — > 00, to be used later. At 
leading order in J Eq. (^) can be reduced to the second order differential equation the solution 
to which is given by 



where?] 



1 



V2 



T 



1 + 21AW 



Wj{v 3l V 2 ) 



1/2 



P 



(^3,^2) 



(2r 2 



(39) 



4 Three-site open spin chain and integrability. 

As we have established in the preceding sections the original problem of calculation the spectrum of 
the evolution equation (|6]) is reduced in multicolour limit to the diagonalization of the Hamiltonian 
( |i8|) which defines the one-dimensional open spin chain with three sites. Moreover, the operators 

~2 *• ~2 

J and J 3 are shown to be the integrals of motion of the system, i.e. [(«/ , J 3 ),7i]_ = 0. In 

9 Here and below the powers of factor of J in front of functions in continuous limit are introduced in order to 
have appropriate scaling properties for J — > 00 when X}-f=o ~~ * J Jo 
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addition to that there exists another 'hidden' conserved charge j22] which thus leads to an exact 
integrability of the three-particle system: the number of integrals of motion equals the number of 
degrees of freedom. Let us show that the system described by the Hamiltonian is equivalent to 
an integrable one-dimensional open spin chain with impurities. For these purposes we will use the 
Quantum Inverse Scattering Method (QISM) formalism. The latter has been used earlier to prove 
the integrability of the multi-reggeon interaction in QCD relevant for the high-energy behaviour 
of physical cross sections and identify the corresponding model with the XXX Heisenberg magnet 
of non-compact spin s = [[25|, p6[| . 

The central object in the algebraic operator formalism of the QISM |28| is the i?-matrix 
which depends on a complex spectral parameter A and acts on the tensor product of vector spaces 
V a ®V b having the dimensions of local spaces on each site. The indices a and b refer to the quantum 
and auxiliary spaces, respectively. The i?-matrix obeys the star-triangle Yang-Baxter equation 

mm 

R a , b (X - fi)R a>c (X)R b M = R b MR a , c (X)R a>b (X - /i). (40) 

In our consequent considerations we only discuss the three-site open spin chain. The mon- 
odromy matrix for it is defined as usual by 

T b (X) = R ai ,b(X - 5i)R a2:b (X - 5 2 )R a3>b (X - 5 3 ), (41) 

where the impurities, Si, are introduced in order to enrich the family of representations. The 
transfer matrix for open spin chain ]3D| 

4(A) = tr b K + (X)T b (X)K-(X)T b -\-X), (42) 

with trace taken over the auxiliary space, — has been shown by Sklyanin to form a family of 
mutually commuting operators [^9| 

[4 1 (A 1 ),4 2 (A 2 )]- = 0. (43) 



This is a direct consequence of the Yang-Baxter equation. The boundary matrices in Eq. (|42|) sat- 
isfy the reflection Sklyanin equations |29| and their simplest solution which leads to the conformal 



invariant transfer matrix is given by the unit matrix^ K ± = i. 

For the two-dimensional auxiliary space (corresponding to spin |) the i?-matrix defines the 
Lax operator 

1 x / \ I 73 T- \ 

(44) 




10 Note that similar choice was made in the context of open spin chains arisen in the problem of integrability of 
quark-gluon reggeon equations |30| . 
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with = =F (J 1 ± iJ 2 ^j. The auxiliary transfer matrix which generates the set of local integrals 
of motion can be redefined according to |29[ 



as 



i) = nA- 1 K,i(A-*-l)}«i(A), (45) 

where A {i? a ^i (A)} = A(A + 1) — J 2 stands for the quantum determinant ^2) and 
ti (A) = tri L ai (X - 5 1 )L a2 (X - 6 2 )L a3 (X - S 3 )a 2 Ll 3 (-X - <5s)2£,(-A - 6 2 )L\ i {-X - 6 x )a 2 , (46) 



which is an even function of the spectral parameter [^9| ti(— A) = ti(A). The expansion of this 
expression w.r.t. the rapidity, A, defines the set of commuting integrals of motion. 

For our purposes it is enough to restrict to the following values of the shifts: 61 = 63 = 6, 
62 = 0. An explicit calculation leads to the result 



ti (A) = J?(A) - 4A 2 (A 2 - jl) f - 2X 2 Q, 
where i?(A) is a c-number function 

0(A) = -2 (A 2 + jl - 5 2 ) (A 2 + J 2 ) (A 2 + 51 ~ 6 2 ) + 4A 2 (A 2 - 5 2 ) f j 2 + 4A 2 (fjl - 6 2 g 



(47) 



i2 



J is the total Casimir operator (|25|) and Q is a non-trivial integral of motion 

Q = [J 2 12 , J 2 23 ]+ - 26 2 { f l2 + J23 



(4J 



The fundamental Hamiltonians which commutes with the above set of charges can be con- 
structed following Sklyanin JHJ. The fundamental monodromy matrix is defined by product (|4l|) 
of i?-matrices with equivalent representations of conformal spin in quantum and auxiliary spaces 
and the shifts of spectral parameters defined as above. Then the Hamiltonian of the chain is 
defined as 



n 



1 d 

2 dX 



In 4(A). 



(49) 



A=0 



Due to the diagonal form of the boundary matrices the latter is given, up to irrelevant c-number 
term, by the sum of two-site Hamiltonians |2{| 



7~t = h aia2 + h 



(50) 



which are expressed by the equation 



h. 



ab 



d_ 

dX 



\nR a , b (X-6)=R-l(-6)R'(-6). 



(51) 



A=0 



13 



Making use of the explicit form of the Yang-Baxter bundle p3|, |34 



i^(A) = P ab l { /; b + % (52) 

one easily obtains 

^ = HJab -5)+ ?P(j ab + 5)- 2^(1), (53) 

where we have legitimately added a c- number term to mimic the QCD pair- wise Hamiltonian flT8|). 

The commutativity of Q (]48[) and 7^ (|50|) can be checked a posteriori by an explicit calculation 
of the commutators [Q, ^(1,3)2]-, e.g. in the conformal basis (|^). Here the problem is reduced to 
the evaluation of the only matrix element which reads 

[Q, h 12 ]j,j+l = W 12 • ,j • I + [Jl2\j,j - 2 <5 2 j {[hl2\j+l - [hn]]} [J^L'J+l 

{•^2 ~2 1 "2 

[«^12]j+lJ+l ~~ ["^Lj j [«^23]jJ+l> 

where in the last line we have used the explicit form of the matrix elements of the Hamiltonian 
(|53|). This result together with similar calculation for [Q, ^23]- in Pi3-transformed basis gives 

[Q, = 2[j1 3 , J* 2 ]_, [Q, h 23 ]_ = 2[J* 2 , J 23 ]_, (54) 

and thus [Q,ft]_ = 0. 

Finally, if we will put 5 = | as suggested by QCD (|i~8|) we immediately obtain the conserved 



charge 'empirically' found in Ref. [22 . 



5 Master equation. 

Thus as has been said above the existence of the additional charge Q leads to complete integrability 

of the quark-gluon system (z/ x = z/ 3 = 1, z/ 2 = 2). This implies that the Hamiltonian 7i is a function 

-2 

of the integrals of motion only, J and <2, and allows to reduce the complicated eigenvalue problem 
for the Hamiltonian to the more simple one for the charge Qt 

Q T # = qT $ ^th Q T = [f 12 , f 23 ] + - § {-^12 + ^23} , (55) 

making use of the conformal invariance of the system. The latter condition means that it is 



convenient to look for the solution to (55) in the form of expansion w.r.t. the three-point basis 



H) elaborated in the previous sections, i.e. 



* = X>iPjtf- (56) 

3=0 
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~2 

The matrix elements of the charge possess, similar to the non-diagonal Casimir operators Jn 2 ) 3 , 
only three non-zero diagonals 

{V^A^MQtIVj-A^W) = 6 J> J (WQtIjj + 5 jW [Qr) j+ltj + 5 fjj \ Q-i ,,) , (57) 
where 

[Qt]jj = 2 [fuhAh ~ I {[Jnhj + [Jlsh] , (58) 

[QT) j+1 j = {[JnUw + [Jnh - 1} (59) 

and [Qt]j + i,- = [Qt] jj+ i- Therefore, the above equation ( [55}) leads to a three-term recursion 
relation which can written in a concise matrix form by introducing the two-dimensional vector 

= Mj&j with Mj = f [ St W+i j _ ( 60 ) 

It must be supplied with the boundary conditions = l^j+i = which manifest the polynomial 
character of the eigenfunctions. The solution to the Eq. (|60"D can be presented thus in the form 

= MjMj-i . . . M #o, 

with the 'vacuum' state \Pq — Cq) where we have voluntary assumed $b = 1 due to its arbitrary 
normalization. The boundary condition = presumes that the qx values are the real roots 
of (yet unknown) orthogonal polynomials. For the expansion coefficients themselves it implies 
that the eigenvalues of the matrix Mj must be complex because otherwise the solution will be a 
monotonic function of j on the whole interval j G [0, J] and the only way to fulfill the boundary 
condition is to accept necessarily the trivial solution \Pj = 0. Therefore, from the eigenvalue 
equation, det[Af — M ■ 1] = 0, we get the condition for existence of the non-trivial solution to the 
recursion relation 

([QtIjj ~ qrf < 4[Qt], j+1 [Qt] jj _ 1 . (61) 
Let us redefine the expansion coefficient according to 

*i = (62) 

where 



(j + 1) 3 (j + 3) 3 /t ., 1 w / ,.,^' 

^-—y 3 (J- J + l)(J + J + 5 ) 



(63) 
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Figure 2: The spectrum of the conserved charge qr and a few trajectories from different sets (see 
the text). 

The recursion relation then takes more transparent form which reads 



This master equation will be the main object of our consequent analysis. Unfortunately, there is 
little hope to solve Eq. fl64]) analytically for arbitrary values of the total conformal spin, therefore, 
below we develop a WKB-type expansion for the conserved charge and the energy of the three-site 
chain at large J. 

6 Quantization of integral of motion. 

From the constraint ( |B~T1) one can readily derive the critical values of the spectrum of the conserved 
charge for asymptotical values of the total conformal spin J. This leads to the inequality 



and to the fact that the maximum of the spectrum is achieved for j max = -j^J ■ This means that in 
the vicinity of the upper boundary of the (^-eigenvalues the coefficients #j are peaked around this 



(2j + 3)7}+! + (2j + 5)7}^ + q) (2j + 3)(2j + 5) ([Q T ] jd - q T ) 7} = 0. 



(64) 



< q T /J 4 < ~ 



(65) 
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point. The validity of the above estimate can be checked by explicit numerical diagonalization 
of the matrix (|57f) which gives the quantized values of the charge Qt- The resulting spectrum 
is shown in Fig. |2|. We can immediately notice that the values of qx form the families of one- 
parametric curves. Obviously, there are two possible sets of trajectories which apparently lead to 
an equivalent description of the spectrum. The first set which starts from the top of the spectrum 
(shown by dashed lines in Fig. |2|) and behaves as 

= J* (4°P( n ) + c fP(n) J- 1 + 0(J~ 2 )) . 

These trajectories are labeled by an integer n. Comparing this expression with Eq. (|65|) we 
conclude that CQ 0p (n) = |. Another alternative is to adhere to the parametrization of the curves 
starting from the bottom (shown by solid curves) 

labeled by m. 

These results can be viewed as different analytical continuations of the charge, defined for 
integer J, to complex values of the conformal spin. Whatever choice we made, — that can be 
judged by employing external physical arguments, — the result is the same: the fine structure of 
the spectrum is generated by the pre-asymptotic corrections in the conformal spin J. 

6.1 Description from below. 

Let us address first the question of the description of the spectrum of the charge starting from 
below. In this case the trajectories behave as qx oc J 2 and one can immediately find two exact 
solutions to the master equation 

53 53 
#ac t-i ( J} = — _ + (J + I) 2 , qr ct - 2 (J) = ~ + ( J + 5) 2 , (66) 

separated by a gap from the rest of the spectrum. At this boundary of g^-values the classically 
allowed region, as deduced from Eq. (|6^), tends to spread over the whole interval of attained 
values of j with exception of small vicinities of the end points j en d ~ 1, J. 

To find the WKB solution let us introduce the function Tj = J(— l) J 0(r) with j = tJ. Then 
the recursion relation can be reduced in leading order in J to the second order differential equation 
for 0(r): 

r 2 (l - r 2 )0"(r) - r{\ - r 2 )0'(r) + 2(g* - 1 + 4r 2 )0(r) = 0, (67) 
where q T = J~ 2 qx- The solution to this equation is given by 

0(r) = C^+V+^^^tIt 2 ) + ^-V-^^H^Ir 2 ), (68) 
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with being complex constants and 2 ^ : 'i(±^t|t 2 ) 



21 1 V l±2i77 T 



r 2 ), where we have 



introduced the shorthand notation t]t = ^V^Qt ~ 3. 

To determine the unknown coefficients one has to find the solution of the master equation 
in the vicinity of the deflection points r — > 0, 1, i.e. J 3> j ~ 1 and J> J - j ~ 1, respectively. 

For J ^> J — j ~ 1 the master equation simplifies into 



T j+1 + Tj^ + 2T 3 = 0, 
with boundary condition Tj + i = 0, and the solution to which is 

r, = (-i)V-j + i). 



(69) 



(70) 



An important property of the solutions we have found above is that there exists a region, J — j ~ 
\fj \ where both of them are valid so that by sewing them together we get the unknown coefficients 
CW. Namely, 

C {+) _ 2 .Ti(-ttfr|l) m] 
C(-) ~ *Fi(«*r|l) ' 1 J 

The quantization condition for the charge qx can be obtained from the similar procedure used 

at the other end of values attained by j. However, instead of an explicit solution of the resulting 

reduced (but still quite complicated) recursion relation one can adopt the strategy based on the 

definite symmetry properties of the system under permutation of quark fields. Assuming the 

eigenfunction of the 'hidden' charge to be simultaneously the eigenfunctions of the permutation 

operator P 13 : P 13 V = e^V with up = 0, tt and Pi 3 Pj y (#i, 2 \9 3 ) = V J;j (9 3 , 2 |0i); we get 



k=0 



(72) 



with the phase being the function of the conserved charge qr the dependence on which is deter- 
mined by the relation 



V^ = arg(-l) J E(-i: 



j=0 



(j + l)(j + 3) Tr 



(73) 



Therefore, we can find the boundary solutions making use of known WKB ones. Thus for the 
region J ^> j ~ 1 we get immediately 



r,- 



(j + l)(j + 3) 



U + 2) 

After simple algebra it can be cast to the form 



3/2 



1 dr _ w i( r )0( r ) 



o [r 3 (l -r 2 )]V2 



(74) 



T = ( _ 1)J+J (i + l)(i + 3) ^ 
3 1 J (J + 2) 



x 



C7(+) 



-l) J 2^i(^ T |l) + 



r(| + ^ T )r(|+j-^ T 
r(|-^ T )r(|+j + i?7 T ) 



4^3 (^T 1 1) 



+ + 



• (75) 
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where 4-7-3 (z?7t|1) = 4F3 f ^3^ , +i ^V^^^ lT -l_^^ r lj an d possesses definite scaling properties 
when expressed by the series. Moreover, Eq. ([75]) obviously fulfills the boundary condition T_i = 0. 
From Eq. (|6"8"D it follows that for large J 

T,j = {(7W j 1+2ir > T J~ 2ir > T + cHji-^wj^J, {1 + 0(j/J)} , (76) 
while from ( |75| ) we have 

j\ = ( 1) J+j c * y [c t(+) jf 1 ~ 2i7 '' r r ^ + + ^(-)ji+2ir? T r (2 _ %r l T 



r (I - Z77 T ) ' r (I + 

+ [C (+) 2 ^i(^t|1) + C(-) 2 ^(-^ T |i)] I {1 + 0(1/ j)} . (77) 

Matching these equations in the region of their overlap j ~ \fj making use of Eq. ([H]) results 
into the following quantization condition for the charge 



2?7t m J — arg — ^ r- — — — — — - = 7rm, (7c 



r 


[1 + ivr) 




1) 


r 


[1 ~ iVr) 


2 ^i(i?7r|l) 



with m G Z + . The trajectories generated by this equation are the ones shown by solid lines in 
Fig- 0- 

6.2 Description from above. 

For the description from above the classically allowed region is concentrated in the vicinity of the 
point jmax = 775 Therefore, we expand the master equation around it, j • = 4g (j + Xy/J^j, and 
look for the solution to Eq. fl6"3]) in the form of series w.r.t. inverse powers of J 

00 00 
*(A) = £%)(A)J-* /2 , and g T (J,n) = J 4 £ g^(n) J^, (79) 
£=0 1=0 

where Tj = <£(A). This leads to the following sequence of differential equations 

Z>(i)* ( o)(A)=0, (80) 
V {1) <P {1) (\)+V {2) <P m (\) = 0, (81) 
P(i)^( 2 )(A) + X>(2)0(i) (A) + D(3)0(o)(A) = 0, 

D(1)0(3)(A) + P( 2 )0(2)(A) + Z>(3)#(1)(A) + P(4)^(0)(A) = 0, 
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supplied with boundary conditions $(±00) = and explicit form of differential operators given 
by 

d2 ■ r- ,(i) ox2 



2> (1) = — + 4(6-g^-2A' 
V {2) = ~^ + 8A (6-4V2-A 1 

V {z) = + 4 (3 - A 2 ) + A-^ - (43 - 96^ + 4g? - 24A 2 + 48V2A 2 + 2A< 

o u A ^ ' dX dX ^ 

P (4 ) = -~ + 4A (6 - 4v/2 - A 2 ) ^ + (-12 + 2^ + 3A 2 ) ^ - 2A (57 - 48^ + 8^2A' 
The solution to Eq. (^) is expressed by Hermite polynomials 



$ (0) (A) = # n (Vav^Aj e-^ A2 . (82) 

The solutions to others ([81]) are given by 

* W (A) = / dX'G(X, X')J {t) {\') + C (e) $(p)(\), (83) 

with the Green function of the homogeneous equation (|8~0"D , "D(i)G(A, A') = 5(X — A'), and the 
source = — Xl=o £fy+i-fc)^(fc)- Explicitly they are given by a linear combination of the 
Hermite polynomials, e.g. for the first non-trivial solution we have 



2 



#(i) (A) = — _ £ c fc # fc V 2 V2X , (84) 
2^2^ fc v 7 



with non-zero coefficients c n+3 = — ^, c n+ i = —^(33 — 24v^2 + 3n), c„_i = n(17 — 12a/2 + |n) 
and c„_3 = |n(n — l)(n — 2). Expressions for higher corrections are more involved. 

The quantization for the charge q? immediately follows from boundary conditions. The 
quantization of higher WKB corrections, q^ (£ > 2), appears from the condition of absence of 
zero-mode in the source terms. In this way we derive a few next-to-leading corrections to the 
charge 



1 

2' 



q^\n) = 6--^=-V2n, 

q^\n) = ^1(371 - 72^) + - (ll -72V2)n + y™ 2 , (85) 



(with n — 0, 1, . . .) which allow to describe the upper part of the spectrum with a high accuracy. 
The integer n has an obvious interpretation as the number of nodes of the solutions in the classically 
allowed region. 
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7 Energy. 



The main advantage of having the integrability of the system is that the eigenf unctions of inte- 
grals of motion simultaneously are the ones for the Hamiltonian. By taking the average of the 
Hamiltonian w.r.t. its eigenfunctions ^ we readily obtain 

U{J,n) = 2(j: e )A Ee(j)g^- 3 -, (86) 

with e(j) = + 1) + + 4) + 2 lE . 

At the same time using the representation of the fundamental basis in terms of the hyperge- 
ometric function (|3~3"D and taking into account definite symmetry properties of the system under 
permutation of quarks one can derive the equation linear in the expansion coefficients 

Both of these formulae display the complementary aspects of the formalism. When supplied by 
the solution of the master equation they provide the exact solution of the three-body problem 
under study. 

Thus in order to find the explicit form of the energy as a function of the conserved charges one 
has to substitute the WKB solutions found before into Eq. (^) or (^). 



7.1 Description from below. 

The lowest trajectories corresponding to the solutions (^) were known before || 35, 11 

£ ™,- 1(J) = |{ w+3) + 27E _i__l_j, 

£°**- a (J) = y [m J + 3) + 2 7 e - \ + -j^} . (88) 

The first (second) energy corresponds to the anomalous dimension of the genuine twist-3 part 
of the structure function e {h£) in multicolour QCD since as has been found in the references 
alluded to above the eigenfunctions corresponding to Eq. fl88|) coincide with the coefficient function 
which enter into the relation between e and Hl and corresponding three-particle correlators (@|D. 
Therefore, the description from below seems to have more physical grounds to be trusted as a 
true analytical continuation of the anomalous dimensions to complex J since otherwise it looks 
dubious that the trajectories (|88|) pick out just one point for a given J from another set of 'genuine' 



analytical functions discussed below in the next section |7T2. 



21 



The remainder of the spectrum can be found by inserting the expression for the eigenfunctions 
fl68|) with (|7ID into Eq. flS7|) . We obtain the following result for the energy levels after the gap (cf. 

Jr£(J, rfr) = 2 In J + A lE + 2 Re ^ (J + - | (89) 

As will be shown by numerical analysis this result is valid with high accuracy for the whole 
spectrum. 



7.2 Description from above. 

To evaluate the energy for the top of the spectrum it is more convenient to use the definition 
The energy can be evaluated consistently in WKB approximation as expansion in inverse powers 
of the conformal spin J as 

S(J,n) = ^^\j) + ±S^(n)J-*y (90) 

Substituting the solutions ([82]), Q84|) into (|86|) and equating the coefficients in front of J~ l we get 
the leading and non-leading corrections to the energy 

£ (0) (J) = 41nJ + 4 7jB -21n2- ^, (91) 
8 {1 \n) = 12 -V2- 2V2n, (92) 



Thus at large J the asymptotic behaviour of the energy of the three-site open spin chain is given 
by 

^-S(J, Or) = In fer/2) + A lE - 1 + 0(J~ 2 ). (93) 
And again as before in Eq. ( |89D the energy is a function of the conserved charges of the problem. 

7.3 Comparison with numerical analysis. 

Let us confront the analytical results derived so far with explicit numerical diagonalization of 
the QCD anomalous dimension matrix. This can be done either by using the expression of the 
Hamiltonian in the conformal basis elaborated before 

J-^ = - |) 5 jk + ^W ej (J)e(£)W ek (J), (94) 

or making use of known expression for the multicolour anomalous dimension matrix in the basis 
of ordinary local operators Zj-j oc ijj(id + y gG + ±{id + ) J ~^il) which is more suitable for numerical 
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Figure 3: The energy spectrum and the trajectories (in dashed) defined by Eq. (p9|) . The lowest 
two eigenvalues (solid curves) are given by (|S8|). 
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Figure 4: The spectrum of energy and a few (dashed) trajectories described by the analytical 
formulae of section |7.2| . 
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handling p5L ITT 



w H jk - 5 jk (e(j) - - A ) - 9 ^ {j - k){j + 2 )( j + 3) + {k^J-k)- (95) 

Here we have used the following step functions = {1, if j > k] 0, if j < k}. 

The result of numerical diagonalization is presented in Figs. ^| and [| together with trajectories 
derived from the study of the master equation. As we have mentioned above the Eq. fl89|) gives a 
perfect description for the entire spectrum of the anomalous dimensions. While for the top of the 
spectrum, due to limited number of terms taken in the WKB expansion of the energy, we cannot 
penetrate too far inside the spectrum since otherwise the condition J>n, — whose validity was 
assumed in the reduction of the recursion relation to the set of differential equations (|80| , |8T|) , - 
will be violated. 



8 Conclusions. 

Making use of the algebraic operator formalism of the Quantum Inverse Scattering Method we 
have identified an integrable one-dimensional inhohmogeneous three-site open quantum spin chain 
model to describe the evolution of the twist-three chiral-odd quark-gluon correlators in Quantum 
Chromodynamics in multicolour limit. The integrability of the evolution equations results into 
substitution of the original complicated problem of diagonalization of anomalous dimension matrix 
to the more simple one for the 'hidden' conserved charge. However, since the problem turns 
out to be still quite complicated to derive the analytical solution we have employed the WKB- 
type expansion w.r.t. the total conformal spin of the system. Taking the first few non-leading 
corrections allows to describe the energy spectrum analytically with good accuracy which can be 
systematically improved, however, at a price of a rather involved analysis. The 'hidden' charge 
unravel the structure of the spectrum — it distinguishes the particular components of three-parton 
correlators with different renormalization group dependence. 

Once the 1/N% corrections are taken into account the integrability of the system will be im- 
mediately violated. An interesting question is how the trajectories of multicolour anomalous 
dimensions will mix with each other once the 1/N% effects are switched on. 

The same techniques, up to minor modifications required due to the loss of definite permutation 
symmetry of quark fields, can be applied to the chiral-even sector corresponding the transverse 
spin structure function gi{x). 
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